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ABSTRACT 

He  develop  a  mathematical  model  of  the  ear's  response  to  Intense  sounds 
(over  120  dB)  that  starts  with  pressure  In  the  free- field,  aocounts  for  the 
effects  of  the  head,  external  ear,  middle  ear,  and  ooohlea,  and  ends  with 
predictions  of  basilar  membrane  motion.  The  model's  structure  parallels  the 
ear's  anatomy  and  funotlon  and  fits  hearing  loss  data  for  weapons  Impulses 
that  are  not  explainable  by  any  other  meohanlsm.  The  model  suggests  that 
upward  movements  of  the  basilar  membrane  produoe  the  stresses  that  are  the 
primary  cause  of  damage. 


IETRODUCTIOH 

The  previous  paper  In  this  symposium  pointed  out  deficiencies  In  present 
damage-risk  criteria  (DRCs)  whloh  ranged  frcm  a  lack  of  theoretical  foundation 
to  outright  inaoouraoy,  especially  where  large  oaliber  weapons  are  concerned 
(Prloe,  1986b).  This  analysis  was  based  In  part  on  a  series  of  experiments 
suggesting  that  the  ear's  response  to  intense  stimulation  oan  be  understood 
best  as  a  function  of  mechanical  stress  within  the  Inner  ear  (Price,  1979; 
1981;  1982;  1983a;  1983b;  19830;  1986a).  However,  we  can't  measure  stress 
within  the  inner  ear  dlreotly;  therefore  mathematical  modeling  of  the  ear's 
behavior  at  high  intensities  is  the  most  viable  alternative. 

Such  a  model  is  the  subjeot  of  this  paper.  Given  that  the  participants 
in  this  symposium  have  technical  backgrounds  In  areas  other  than  cochlear 
modeling,  we  will  try  to  present  enough  detail  to  make  the  Implications  of  the 
model  Intelligible.  For  those  interested  In  the  complete  development,  we 
intend  to  publish  a  fuller  aooount  In  an  appropriate  publication.  In 
developing  this  model  there  were  a  number  of  general  considerations  that 
guided  the  modeling.  First,  there  was  the  desire  to  be  able  to  relate  the 
various  elements  of  the  model  to  specifio  physiological  entities  and  to  give 
them  physically  reallstlo  values.  This  procedure  has  heuristic  value  and  has 
a  maximum  chance  of  providing  physical  insight.  Furthermore,  this  procedure 
should  allow  modeling  of  the  ears  of  more  than  one  species,  eaoh  with  Its  own 
characteristic  structures  and  values.  Such  a  capaolty  is  desirable  because 
certain  kinds  of  data  are  available  on  only  some  species  because  of 
anatomical,  practical  or  ethioal  considerations.  Secondly,  the  model  should 
be  calculable  with  modest  oomputer  resources.  Third,  the  main  foous  of  the 
model  was  on  the  behavior  of  the  ear  at  high  Intensities  where  loss  was 
thought  to  be  a  funotlon  of  meobanlcal  stress  C 1 30  dB  and  up).  And  lastly, 
the  model  had  the  goal  of  taking  free-fleld  sound  pressures  as  the  aooustlc 
input  and  passing  the  energy  through  the  outer,  middle  and  Inner  ears  to 
ultimately  provide  an  estimate  of  stress  within  the  sensory  oells  of  the  Inner 
ear  (the  organ  of  Corti).  These  last  goals  required  that  models  of  the  head, 
outer  ear,  middle  ear,  and  Inner  ear  all  be  integrated  and  allowances  made  for 
non-linear  elements  In  them. 


137 


Assuming  that  the  foregoing  efforts  would  meet  with  some  suooess,  we 
hoped  to  'explain*  certain  perplexing  data  on  the  ear's  susceptibility. 
Specifically,  researoh  In  two  different  laboratories  has  shown  the  ear  to  be 
most  susceptible  to  Impulses  with  speotral  peaks  In  the  mid-range  and  muoh 
less  susoeptlble  to  Impulses  with  their  peak  energy  at  lower  frequencies 
(Dancer,  Lenoir,  Buck  and  Vassout,  1983;  Prloe  1983a;  1986b).  In  experiments 
with  the  oat  ear  and  Impulses  from  weapons,  howitzer  Impulses  were  found  to  be 
as  hazardous  as  rifle  Impulses  only  when  the  howitzer  Impulses  were  about  10dB 
more  Intense  (Price,  1983a;  1986b).  Furthermore,  losses  tended  to  be  In  the 
middle  range  of  frequencies  for  both  types  of  Impulse  (Prloe,  1983a;  1986b; 
Price  and  Llm,  1984).  These  findings  are  cruolal  to  the  rating  of  hazard, 
beoause  all  DRCs  would  have  predloted  essentially  the  opposite  result. 
Furthermore,  no  traditional  measures  of  sound  (energy,  A-weighted  energy, 
eto.)  produoe  the  same  ordering  of  effeot. 

In  this  paper  we  will  describe  the  models  used  for  each  of  the  elements 
In  this  Integrated  chain,  we  will  compare  the  model's  output  with  appropriate 
data,  and  lastly,  we  will  compare  the  model's  output  with  various  weapons 
Impulses. 


THE  MODEL 

In  describing  the  model,  the  details  presented  will  be  appropriate  for 
the  cat  ear.  This  choice  is  made  beoause  the  physlologioal  and  aooustlc 
values  are  best  known  for  the  oat  and  also  because  the  hearing  loss  data  to  be 
explained  were  produoed  with  the  oat  ear.  The  ears  of  other  mammals  are 
similar  and  the  principles  used  In  the  modeling  should  transfer  with  only 
modest  adaptations  to  fit  anatomic  details.  A  schematic  diagram  of  the  human 
ear  with  the  anatomical  elements  labeled  Is  provided  In  Fig.  1  as  an  aid  in 
following  the  development  of  the  model. 

Given  the  large  number  of  authors  who  over  the  years  have  proposed  models 
for  various  parts  of  the  auditory  system,  it  is  difficult  to  parcel  out  the 


Vestibular  volume 


Fig.  1 .  Schematic  diagram  of  the  ear. 
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intellectual  debt  accurately.  In  producing  this  synthesis  the  authors 
gratefully  acknowledge  the  efforts  of  the  many  who  have  preceded  them  and 
apologize  for  any  oitations  that  are  missed. 

Diffraction  fllr  Vestibular 


Fig.  2.  Schematic  diagram  of  the  ear  model.  P  is  the  free-field  sound 
pressure;  Pe  is  the  pressure  at  the  eardrum;  Ue  is  the  volume  velocity  into 
the  eardrum.  Pc  is  the  pressure  aoross  the  BM  at  the  base;  Uc  is  the  volume 
velocity  into  the  cochlea  at  the  oval  window.  The  values  of  the  elements  for 
Bulla  closed  are:  Rdf  =  rho*c/Shead;  Ldf  =  0 .5*rho/SQRT(PI*Shead) ;  rho  = 
1.13E-3  gm/cnT3;  c  =  3.52E4  om/sec;  Shead  =  PI*Rhead*Rhead;  Rhead  =  1.65  cm; 
Rpl  =  rho*c/S1 ;  Lpl  =  0 ,85*rho/SQRT(PI*S1 ) ;  LI  s  1.185  cm;  L2  =  L3/3;  L3  = 
1  .705  cm;  SI  *  5.065  cmA2;  S2  =  2*S3;  S3  =  0.375  cm“2;  Lds  =  0.0482;  Cds  = 
4.70E-7;  Rds  =  715;  Ldm  =  0.0214;  Cdc  =  3.62E-7;  Hdc  =  60;  Rh  =  20;  Lh  = 
0.01348;  Cb  = '3.806E-7;  Ct  s  1  .875E-7;  Nt  =  53 r  Cmi  =  7E-12;  Rmi  =  6.4E5;  Li  = 
1.6;  Cis  s  8.4E-12;  Rls  =  1E4;  Ls  =  3.3;  Cal  s  1 .1  IE-9;  Hal  =  2E5;  Lv  =  22;  Ro 
s  1.2E6;  Ro  =  2.8E5;  Lo  =  2250;  Crw  =  1  .0E-8;  PI  =  3.1415.  In  case  the  bulla 
ls  opened  but  the  bony  septum  is  not  removed,  three  elements  change  in  the 
bulla  model:  Cb(open)  s  Cb#4.70;  Ct(open)  =  Ct*1.06;  and  Lh(open)  s  Lh*1.l6. 
In  case  the  bony  septum  is  removed  also,  then  the  bulla  Impedance  beoomes 
negligible  and  is  short-circuited  in  the  model.  The  units  of  all  masses 
(induotors)  are  gm/om~4;  viscous  loss  elements  (resistors)  are  dyne-sec/om*5 ; 
and  compliances  (capacitors)  are  cm*5/dyne.  The  formulas  for  the  tube 
sections  in  the  outer  ear  are: 

PQ  =  PjOos(kL)  +  i  0iZosin(kL) 

U0  =  Uicos(kL)  +  i  PjZc-^Bin(kL) 


where : 


Pi  =  pressure  magnitude  at  tube  input  end. 

s  volume  veloolty  magnitude  at  tube  input  end. 

P  s  pressure  magnitude  at  tube  output  end. 

0  =  volume  veloolty  magnitude  at  tube  output  end. 

Z°  =  rhoBo/S  ,  characteristic  impedanoe  of  tube. 

S  =  cross  seotional  area  of  tube. 

L  =  length  of  tube. 

k  *  2*PI*f/c  ,  wavenumber  of  waves  in  tube, 
f  =  frequency  in  Hz. 
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Fig.  2  shows  both  a  sohematio  diagram  of  the  ear  model  along  with 
acoustic  and  electrical  quantities  whloh  can  be  measured  at  each  location.  At 
each  stage  the  energy  is  filtered  and  transformed  as  it  passes  through. 
Acoustic  vibration  in  the  outer  ear  becomes  meohanioal  vibration  in  the  middle 
ear,  whloh  becomes  fluid  vibration  in  the  inner  ear,  whloh  ultimately  becomes 
eleotro-chemical  energy  in  the  sensory  cells  and  nerve  fibers. 

From  The  Free-Field  to  the  Eardrum 

The  head  is  modeled  as  a  spherical  baffle  around  which  the  sound  wave 
diffracts  and  in  whloh  the  entranoe  of  the  earoanal  is  located.  The  field 
surrounding  the  head  and  in  the  concha  (the  major  cavity  in  the  external  ear) 
is  approximated  by  a  simple  network  proposed  by  Bauer  (1967).  Next  in 
sequenoe  is  Wiener,  Pfeiffer  and  Baokus*  (1965)  two  tube  model  of  the  earoanal 
which  carries  the  energy  as  far  as  the  eardrum.  Component  values  (Fig.  2)  in 
the  model  were  adjusted  to  fit  the  data  from  Wiener  et  al.  (1965)  on  the 
pressure  transfer  from  the  free- field  to  the  eardrum  (Fig.  3a).  The  agreement 
between  the  model  and  the  data  is  good  for  frequencies  up  to  10  kHz.  We 
believe  the  lack  of  agreement  at  higher  frequencies  is  due  to  measurement 
difficulties  Inherent  at  high  frequencies. 


Fig.  3*  Comparison  of  model-oaloulated  values  with  measurements  of  four 
acoustical  quantities. 
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i. 


Across  the  Middle  Ear 


The  essential  function  of  the  middle  ear  Is  to  match  the  lev  Impedance  of 
air  to  the  high  impedanoe  of  the  ooohlea.  The  model  at  this  point  Is  based  on 
one  by  Zwislookl  (1962)  with  additional  Impedance  elements  given  by  Lynch 
(1982).  The  eardrum  Is  a  complex  structure  whloh  does  not  move  as  a  piston 
but  actually  vibrates  In  segments.  The  part  of  the  eardrum  that  conducts 
sound  to  the  malleus  Is  represented  by  a  series  arm  while  the  part  which  moves 
Independently  Is  described  by  a  shunt  arm.  Losses  at  the  ligaments  Joining 
the  ossicles  are  also  represented  by  shunt  elements.  A  transformer  models  the 
acoustic  transformations  due  to  lever  action  of  the  ossicles  and  the  ratio  of 
areas  of  the  eardrum  and  stapes.  The  connected  volumes  of  the  bulla  located 
behind  the  eardrum  oontribute  to  Its  stiffness  and  produce  a  resonance  and 
anti-resonance  near  4  kHz  described  by  a  parallel-series  element.  Experiments 
are  often  conducted  with  the  middle  ear  opened.  In  this  case,  Increasing  the 
value  of  the  bulla  elements  (Fig.  2)  allows  the  model  to  aooount  for  the 
middle  ear  change.  The  annular  ligament  which  surrounds  the  stapes  at  the 
entrance  of  the  cochlea  is  described  by  a  vlsoo-elastlo  series  element.  The 
fluid  In  the  vestibular  volume  is  given  by  a  mass  element  while  the  round 
window  is  modeled  as  a  compliance  element.  Finally,  the  fluid  load  of  the 
cochlea  and  the  low-pass  pressure  relief  of  the  hellootrema  are  described  by 
mass  and  resistance  elements. 

Model  predlotions  are  compared  with  several  aooustlcal  measurements  at 
the  eardrum  shown  in  Fig.  3b,  3c,  and  3d.  The  prediction  of  the  ratio  of 
stapes  displacement  to  eardrum  pressure  (Fig.  3b)  Is  seen  to  agree  very  well 
with  the  data  from  Guinan  and  Peake  (1967)  •  They  show  a  uniform  response  up 
to  1  kHz  and  a  12  dB/oot  rolloff  at  higher  frequencies.  In  Fig.  3c  we  see 
that  the  model  (solid  line)  reproduces  the  input  Impedance  at  the  eardrum  as 
measured  by  Moller  (1963)  and  Tonndorf  and  Kbanna  (1967)*  At  frequencies 
above  4.0  kHz  there  are  large  discrepancies  between  measures  of  the  transfer 
funotlon  largely  because  of  the  practical  problems  assooiated  with  making  the 
measurements.  The  wavelength  of  the  sound  at  these  frequencies  beoomes  short 
enough  to  be  a  problem,  the  eardrum  Is  really  an  extended  surface  with  a 
complex  shape.  Furthermore,  at  high  frequencies  it  does  not  vibrate  as  a 
single  unit  but  breaks  up  into  oomplex  modal  patterns.  Finding  a  single  point 
to  make  a  representative  measurement  is  not  really  possible.  It  Is  therefore 
not  surprising  that  experimental  measurements  of  aooustlo  values  for  the 
middle  ear  and  external  ear  might  show  considerable  divergence  at  specific 
points  (even  without  considering  the  problem  of  individual  differences  between 
animals).  Fig.  3d  shows  the  ratio  of  eardrum  pressure  to  intracochlear 
pressure.  The  model  matches  Hedzelnltsky '  s  (1980)  measurements  reasonably 
well. 

The  Cochlea 

The  cochlea  is  modeled  as  a  two- chambered,  fluid -filled,  box  with  rigid 
side  walls.  The  partition  separating  the  chambers  Is  rigid  except  for  a 
tapered  basilar  membrane  (BM)  which  is  narrow  at  the  base  and  widens 
exponentially  towards  the  apex.  Its  meohanioal  properties  also  change  along 
its  length;  the  meohanlcal  stiffness  per  unit  area  decreases  exponentially  and 
the  mass  and  viscous  loss  per  unit  area  increase  exponentially  along  its 
length.  The  mass  inoludes  the  membrane,  the  organ  of  Corti,  and  the  fluid 
which  moves  along  with  it.  In  contrast  to  the  basilar  membrane,  the  oross- 


141 


sectional  area  of  the  ooohlea  decreases  exponentially  from  base  to  apex.  At 
any  given  looatlon  along  the  BM,  we  assume  that  the  pressure  is  constant 
everywhere  within  the  cross  section.  This  Is  equivalent  to  assuming  the 
wavelengths  on  the  BM  are  long  o  can  pared  to  the  height  of  the  chambers. 

For  sinusoidal  stimulation  at  the  stapes,  a  traveling  wave  appears  on  the 
BM.  For  eaob  frequency,  the  BM  displacement  reaohes  a  maximum  at  a 
characteristic  location,  high  frequencies  having  their  maxima  at  the  stiffer, 
basal  end  of  the  cochlea  and  low  frequencies  having  their  maxima  at  the  more 
compliant,  apical  end  of  the  BM.  Given  the  presence  of  a  traveling  wave 
within  the  ear,  we  have  solved  for  the  BM  motion  with  a  mathematical  procedure 
known  as  the  VKB  method  (Zweig,  Lipes  and  Pierce,  1976)  whioh  Is  applicable 
for  systems  with  traveling  waves,  provided  the  wavelength  is  short  compared  to 
the  dimensions  of  the  system.  This  limitation  of  the  VKB  method  is  admittedly 
a  problem  as  the  driving  frequency  gets  lower.  However  the  method  offers  the 
overriding  advantage  that  dosed- form  analytio  solutions  are  obtainable. 

It  was  mentioned  earlier  that  the  hypothesized  hearing-loss  mechanism  at 
high  intensities  was  stress  within  the  structures  in  the  inner  ear.  However, 
given  the  limitations  of  knowledge  about  the  miorcmeohanics  of  the  cochlea,  we 
have  aocepted  BM  displacement  as  a  first  approximation.  The  stapes  to  basilar 
membrane  displacement  ratio  has  been  calculated  for  the  7.4  kHz  'plaoe'  on  the 
basilar  membrane  for  tones  with  SPLs  ranging  from  55  to  105  dB  (Fig.  4).  Also 
plotted  for  oomparison  in  Fig.  4  are  Rhode's  (1980)  data  for  the  displacement 
at  the  equivalent  cochlear  place  in  the  squirrel  monkey  ooohlea.  Two  things 
are  apparent.  First,  the  calculated  displacement  ratios  are  a  reasonable 
approximation  of  the  measured  ratios.  Secondly,  that  at  low  levels,  the 
response  is  more  peaked  than  it  is  for  high  levels  of  stimulation.  The  fit  to 
the  data  was  produced  by  changing  the  damping  factor  used  in  the  calculation 
(damping  values  listed  on  the  figure).  A  more  easily  interpretable  picture  of 
what  happens  at  a  particular  cochlear  looatlon  as  intensity  rises  is  shown  in 
Fig.  5  whioh  inoludes  essentially  the  same  Information  as  Fig.  4,  except  that 
displacements  are  plotted  instead  of  ratios. 

Recent  researoh  has  suggested  that  the  there  is  an  'amplifier*  in  the 
organ  of  Corti.  If  this  amplifier  worked  at  low  sound  intensities  by  reducing 
losses  in  the  basilar  membrane  (smaller  damping  factor,  better  tuning)  then 
the  sensitivity  of  the  ear  would  be  enhanoed  when  it  was  needed.  Conversely, 
at  high  intensities,  the  effect  of  the  oochlear  amplifier  would  be  overwhelmed 
by  the  relatively  large  forces/displaoements  and  the  response  of  the  BM  would 
be  a  funotion  of  the  essential  properties  of  the  membrane/ fluid  system  itself. 
Therefore  in  working  with  noise  exposures  and  this  model,  we  have  used  the 
higher  damping  values  appropriate  to  high  levels. 

The  plot  of  predicted  basilar  membrane  displacement  along  its  length  for 
pure  tone  inputs  at  the  stapes,  is  shown  in  Fig.  6.  The  envelopes  of 
traveling  wave  displacements  in  Fig.  6  matoh  essential  features  of  what  is 
known  about  the  ear's  response.  First,  the  looatlon  of  the  peak  varies  as 
the  logarithm  of  the  frequency,  matching  the  map  produoed  for  the  oat  cochlea 
by  Liberman  (1982)  and  'corrected'  by  the  reduced  damping  noted  in  the 
previous  paragraph.  Second,  the  amplitude  of  the  peak  increases  by  3  dB  for 
each  octave  increase  in  frequency.  This  value  is  close  to  that  observed  by 
von  Bekesy  (1949).  Third,  the  Q  of  the  peaks  increases  with  frequency. 
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Fig.  4.  Comparison  of  model  calculation  of  basilar  membrane  to  stapes 
displacement  ratios  at  the  7.4  kHz  plaoe  for  different  driving  frequencies  and 
Intensities.  The  data  points  are  from  Rhode's  measurements  on  the  squirrel 
monkey  (1980), 


Fig.  5.  Model  calculation  of  the  basilar  membrane  displacement  envelope  at 
the  7*4  kHz  place  and  sound  pressures  from  55  dB  to  over  105  dB. 
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Fig.  6.  Calculated  envelope  of  traveling  wave  displacements  for  different 
frequencies. 

The  Stapes 

Thus  far,  the  components  of  the  model  have  been  linear  elements.  At  high 
levels,  however,  the  annular  ligament  of  the  stapes  introduoes  important  non- 
linearities,  (Price,  1974).  The  ligament  is  rugged  and  is  linear  for  almost 
all  sound  pressures  to  whioh  the  ear  is  likely  to  be  exposed  (Wever  and 
Lavrenoe,  1954).  However,  at  very  high  Intensities  the  ligament  limits  stapes 
motion,  effectively  clipping  the  cochlear  input.  Guinan  and  Peake  (1967) 
published  data  on  stapes  motion  at  high  sound  pressure  levels  which  are 
reproduced  in  Fig.  7  along  with  a  non-linear  model  (solid  line)  fitting  the 
data.  The  ligament  model  predicts  that  stapes  displacements  will  reach  an 
asymptote  at  about  42  mlorons  peak  to  peak  displacement.  The  solid  line  of 
the  model  fits  the  displacement  data  very  well.  The  sound  pressure  required 
to  reach  the  displacement  limit  will  of  oourse  be  a  function  of  the  frequency 
of  the  driving  stimulus  and  has  in  essence  been  calculated  for  pure  tone 
stimulation  (Prloe,  1974).  On  the  other  hand,  the  model  makes  it  possible  to 
do  the  calculation  for  spectrally  complex  stimuli  in  which  we  have  a  practical 
interest.  Wever  also  showed  (pp.  148-151)  the  eardrum  moves  linearly  for 
displacements  of  more  than  100  mlorons;  hence  allowing  for  the  middle  ear 
lever  ratio  of  about  2.5,  the  eardrum  would  not  produce  any  additional 
distortion. 

In  applying  the  clipping  calculation  to  the  model,  we  recognize  that  the 
addition  of  a  non-linear  element  must  to  some  extent  be  reflected  in  the 
responses  of  other  middle  ear  structures.  The  present  calculational  scheme 
makes  no  specif io  allowance  for  suoh  effects.  This  may  not  be  too  serious  a 
limitation,  given  that  the  middle  ear  is  critically  damped  and  does  not  show 
appreciable  energy  storage  and  ringing,  at  least  at  the  levels  for  which 
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Fig.  7 .  Comparison  of  clipping-model  output  for  stapes  displacement  with 
empirically  determined  values  from  Guinan  and  Peake  (1967). 

measurements  have  been  made.  The  data  from  Guinan  and  Peak  (1967)  in  Fig.  7 
are  consistent  with  this  assumption.  Furthermore,  in  modeling  the 
otoaolerotic  ear  (a  condition  in  which  the  stapes  is  immobilized  in  the  oval 
window),  Zwlslockl  (1962)  pointed  out  that  the  effect  of  stapes  immobilization 
on  the  impedanoe  at  the  eardrum  is  negligible,  due  to  the  high  impedance  of 
the  lnoudo-stapedlal  Joint  and  the  low  ooupling  Impedance  of  the  eardrum  at 
high  frequencies. 

In  oaloulating  with  the  present  model,  the  solution  for  BM  motion  is 
found  in  three  steps.  First,  the  linear  solution  is  found  for  stapes  motion, 
ignoring  annular  ligament  non-linearities.  Next,  the  linear  solution  is  then 
non-linearlzed  aooordlng  to  the  empirically  determined  time-domain  dipping 
transfer  funotion  of  Fig.  7.  And  third,  the  modified  motion  is  used  as  input 
to  the  ooohlear  portion  of  the  model  whioh  in  turn  prediots  BM  displacements 
at  speoiflc  points. 


Given  stapes  displacement  as  an  intermediate  point,  the  transfer  function 
from  stapes  displacement  to  BM  displacement  at  a  dlstanoe  x  from  the  stapes 
for  a  driving  frequency  f  is  D(x,f).  The  sign  convention  is  that  an  initial 
displacement  (positive)  of  the  stapes  into  the  soala  vestibuli  causes  a 
downward  (negative)  displacement  of  the  BM  (from  soala  vestibuli  toward  soala 
tympani).  The  calculations!  formula  is: 


D(x,f) 

where: 
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=  1 .26  ram2  ,  area  of  stapes  footplate. 

=  0.08  mm  ,  width  of  BH  at  stapes  end. 

=  4.68  mm  ,  tuning  faotor  of  BM. 

=  0.024,  damping  factor  of  BM  at  stapes  end. 

=  distance  from  stapes  In  mm  along  BM. 
a  3.16  t  number  of  traveling  waves  on  BM  at  one  time, 
a  frequency  of  stapes  osoillatlon  In  kHz. 
a  66.1  kHz  ,  looal  resonanoe  frequency  of  BM  at  stapes  end, 

a  f0e“xy,d  ,  looal  resonance  frequency  of  BM  at  looatlon  x. 

deltaAexy,d  ,  damping  faotor  of  BM  at  location  x. 
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The  empirical  formula  for  converting  the  linear  solution  for  stapes  motion 
Into  an  approximate  non-linear  solution  Is: 


xim 

xni  a  ;  "  7 

1  +  *xlln/xmaxl 

where: 

x^^q  a  linear  stapes  displacement  in  microns;  calculated  from  the  linear 
ear-model  with  no  limitation  by  the  annular  ligament. 
xlln  =  ?e*Alf'  In  the  case  of  the  low  frequenoy  tone  used  by  Guinan  and 
Peake. 

Alf  =  °*1^2&-6  cm“3/dyne,  low-frequency  closed-bulla  transfer  ratio  of 
the  middle  ear  listed  in  Table  II  of  Guinan  and  Peake. 

Pe  =  magnitude  of  sound  pressure  at  eardrum. 

xmfiy  a  20.9  microns  ,  maximum  stapes  displacement  extrapolated  from  pure 
tone  data  of  Guinan  and  Peake  assuming  symmetrical  peak  clipping. 
xnl  =  non-linear  stapes  displacement  in  microns. 

APPLICmOH  OF  THE  MODEL 

Work  in  applying  this  model  has  just  begun;  therefore  many  details  of  Its 
application  remain  to  be  explored.  However,  the  early  results  are  promising. 
The  model  does  provide  plausible  explanations  for  the  two  major  'anomalies'  in 
the  hearing  loss  data,  namely,  the  faot  that  for  equal  peak  pressure  levels 
(PPLs)  the  rifle  is  more  hazardous  than  the  howitzer  and  the  fact  that  loss 
occurs  In  the  mid-range  of  frequencies,  even  when  the  spectral  peak  is 
elsewhere. 

The  prediction  of  loss  is  the  mid-range  was  made  some  years  ago  (Prioe, 
1977;  1981)  based  primarily  on  the  tuning  of  the  external  and  middle  ears  and 


146 


the  relatively  flat  spectra  of  commonly  encountered  noises.  Because  the 
present  model  embodies  the  same  transfer  functions)  it  is  also  consistent  with 
the  earlier  contentions.  Later  studies  with  noise  exposures  have  shown  that 
exposures  to  impulses  from  both  the  rifle  and  howitzer  produce  histologically 
and  audiometrically  measured  damage  to  the  middle  of  the  cochlea  (Price) 
1983a;  1986b;  Price  and  Lim,  1984) 

The  more  puzzling  problem  was  the  prediction  of  the  relative  ranking  of 
the  hazard  from  large  and  small  caliber  weapons.  The  explanation  that  has 
resulted  from  work  with  the  model  is  presented  in  the  next  two  figures.  Fig. 
8a  shows  the  stages  in  the  ear's  response  to  a  rifle  at  145  dB  PPL.  This 
impulse  and  the  others  used  here  are  those  actually  recorded  during  impulse 
noise  exposures  and  include  all  the  features  normally  found  in  this  type  of 
data,  such  as  the  ground  reflections.  The  top  panel  shows  the  free-field 
pressure,  typical  of  a  rifle.  The  middle  panel  shows  the  corresponding  stapes 
displacement.  The  largest  displacement  was  to  the  positive  phase  of  the 
impulse  and  was  about  7  microns,  the  ligament  model  clipping  only  slightly. 
The  lowest  panel  shows  the  predicted  BM  displacement  at  the  3  kHz  plaoe  in  the 
ooohlea,  the  area  where  the  largest  effeots  have  been  noted.  Note  that  the 
largest  displacement  of  the  BM  was  about  12  microns  upward  and  is  in  response 
to  the  rarefaotlon  phase  of  the  waveform.  The  parallel  set  of  data  is  shown 
in  Fig.  8b  for  the  rifle  impulse  at  155  dB  PPL.  The  stapes  displacement 
showed  some  clipping  (it  actually  rose  less  than  6  dB  to  13  microns  for  a  10 
dB  lnorease  in  PPL).  On  the  BM,  the  largest  displacement  was  still  upward  and 
was  more  than  25  microns. 


Fig.  8.  Displacements  within  the  ear  calculated  for  a  rifle  impulse  at  145  dB 
and  155  dB  PPL. 
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For  the  rifle,  the  largest  displacements  are  In  the  upward  direction. 
Therefore,  we  might  suspect  that  scala  vestibul 1-ward  movements  do  the  most 
damage.  We  In  faot  suspect  that  hazard  may  well  be  a  function  of  the 
direotlon  of  movement.  Such  a  contention  makes  sense  physiologically.  In 
essenoe,  downward  movement  produces  essentially  compressive  stresses  whereas 
upward  movement  produces  primarily  tensive  deformations.  Tissue  tends  to  fall 
in  tension;  therefore  It  Is  reasonable  to  suspect  that  upward  BM  motions  might 
be  the  most  damaging. 

The  howitzer  Impulse  in  Fig.  9a  is  154  dB  PPL  and  produces  maximum  stapes 
displaoements  of  about  13  microns,  about  the  same  as  the  rifle  at  the  same 
PPL.  On  the  BH,  the  displaoements  in  the  downward  direction  are  about  the 
same  as  for  the  rifle  at  the  same  PPL.  However,  the  upward  displacements  are 
smaller  than  those  for  the  rifle  by  more  than  a  faotor  of  two.  In  Fig.  9b, 
the  oalculated  response  to  the  howitzer  impulse  at  166  dB  can  be  seen  to 
produce  stapes  displacements  of  almost  20  microns  inward  and  a  little  less 
than  15  outward.  If  no  clipping  had  occurred,  the  peak  inward  displacement 
would  have  been  over  160  microns  I  The  clipping  introduced  by  the  annular 
ligament  of  the  stapes  obviously  has  a  major  effect  on  what  gets  transmitted 
to  the  inner  ear  at  high  PPLs.  The  calculated  BM  motions  make  an  interesting 
pattern.  Perhaps  the  most  important  point  is  that  the  maximum  upward 
displacements  are  on  the  order  of  25  or  30  microns,  about  the  same  as  for  the 
rifle  when  it  was  11  dB  less  intense.  Furthermore,  the  largest  displacements 
occur  not  in  response  to  the  fundamental  waveform;  but  to  the  'bash1  riding  on 
it.  It  is  these  smaller  oscillations  that  produoe  large  stapes  movements,  and 
if  their  timing  is  correot,  the  BM  at  the  3  kHz  place  responds.  It  is 
interesting  to  note  that  an  intense  impulse  with  a  long  A-duratlon  would  push 
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Fig.  9.  Displaoements  within  the  ear  calculated  for  a  howitzer  impulse  at  154 
dB  and  166  dB  PPL. 
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Fig.  10.  Upward  BH  displacement  as  a  function  of  PPL  for  howitzer  and  rifle 
Impulses. 

the  stapes  In  and  hold  It  there,  effectively  blooking  transmission  Into  the 
Inner  ear  during  that  period. 

Fig.  10  summarizes  what  we  believe  to  be  the  major  point  of  the  last  two 
figures.  The  upward  displacement  of  the  BM  as  a  funotlon  of  PPL  grows  as  PPL 
rises  for  both  Impulses,  a  little  faster  for  the  howitzer  than  for  the  rifle. 
However,  note  that  for  equal  displacements,  the  PPL  of  the  howitzer  must  be 
higher  than  that  for  the  rifle  by  about  10  dB.  This  is  the  ranking  of  hazard 
that  the  hearing  loss  measurements  have  produoed. 

These  results  are  encouraging  in  that  they  do  seem  to  predict  what  the 
hearing  loss  data  are  showing  and  harmonize  all  the  data.  And  if  we  speculate 
freely,  the  model  could  be  used  as  a  basis  for  a  convenient  'hazard  meter'  for 
any  type  of  intense  sound*  Furthermore,  the  model  has  heuristio  value  in  that 
it  aan  be  used  to  calculate  ways  of  ohanging  the  wave  form  to  produoe  a  safer 
impulse.  What  is  needed  now  is  confirmation  of  the  model's  predictions  in 
real  ears. 

We  intend  next  to  extend  the  model  to  inolude  different  impulses  on  which 
we  have  hearing  loss  data,  to  generate  predictions  of  hazard  for  different 
speoles  to  inolude  man,  and  to  generate  testable  hypotheses  as  a  means  of 
verifying  the  essential  elements  of  the  model.  What  is  reported  here  is  a 
preliminary  effort;  but  if  it  is  supported  by  subsequent  researoh,  it  is  fair 
to  say  that  the  assessment  of  impulse  noise  hazard  will  undergo  dramatlo 
changes. 
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1  DEFENSE  TECHNICAL 
(PDF)  INFORMATION  CTR 
DTIC  OCA 

1  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
IMAL  HRA 

1  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 

1  GOVT  PRINTG  OFC 
(PDF)  AMALHOTRA 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRM  D 
T  DAVIS 

BLDG  5400  RM  C242 

REDSTONE  ARSENAL  AL  35898-7290 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRS  EA  DR  V  J  RICE 
BLDG  4011  RM  217 
1750  GREELEY  RD 

FORT  SAM  HOUSTON  TX  78234-5002 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRM  DG  I  RUBINSTEIN 
BLDG  333 

PICATINNY  ARSENAL  NJ  07806-5000 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  ARMC  FIELD  ELEMENT 

RDRL  HRM  CH  C  BURNS 
THIRD  AVE  BLDG  1467B  RM  336 
FORT  KNOX  KY  40121 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  AWC  FIELD  ELEMENT 

RDRL  HRM  DJ  D  DURBIN 
BLDG  4506  (DCD)  RM  107 
FORT  RUCKER  AL  36362-5000 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRM  CK  J  REINHART 
10125  KINGMAN  RD  BLDG  317 
FORT  BELVOIR  VA  22060-5828 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRM  AY  M  BARNES 
2520  HEALY  AVE 
STE1172  BLDG  51005 
FORT  HUACHUCA  AZ  85613-7069 


1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRM  AP  D  UNGVARSKY 
POPE  HALL  BLDG  470 
BCBL  806  HARRISON  DR 
FORT  LEAVENWORTH  KS  66027-2302 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRM  AR  J  CHEN 
12423  RESEARCH  PKWY 
ORLANDO  FL  32826-3276 

1  ARMY  RSCH  LAB  -  HRED 

(PDF)  HUMAN  SYSTEMS 

INTEGRATION  ENGR 
TACOM  FIELD  ELEMENT 
RDRL  HRM  CU  P  MUNYA 
6501  E  11  MILERD 
MS  284  BLDG  200A 
WARREN  MI  48397-5000 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  FIRES  CTR  OF  EXCELLENCE 
FIELD  ELEMENT 
RDRL  HRM  AF  C  HERNANDEZ 
3040  NW  AUSTIN  RD  RM  221 
FORT  SILL  OK  73503-9043 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRM  AV  W  CULBERTSON 
91012  STATION  AVE 
FORT  HOOD  TX  76544-5073 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  RDRL  HRM  DE  A  MARES 

1733  PLEASONTON  RD  BOX  3 
FORT  BLISS  TX  79916-6816 

8  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  SIMULATION  &  TRAINING 
TECHNOLOGY  CENTER 
RDRL  HRT  COL  G  LAASE 
RDRL  HRT  I  MARTINEZ 
RDRL  HRT  T  R  SOTTILARE 
RDRL  HRT  B  N  FINKELSTEIN 
RDRL  HRT  G  A  RODRIGUEZ 
RDRL  HRT  I  I  HART 
RDRL  HRT  M  C  METEVIER 
RDRL  HRT  S  B  PETTIT 
12423  RESEARCH  PARKWAY 
ORLANDO  FL  32826 

1  ARMY  RSCH  LABORATORY  -  HRED 
(PDF)  HQ  USASOC 

RDRL  HRM  CN  R  SPENCER 
BLDG  E2929  DESERT  STORM  DRIVE 
FORT  BRAGG  NC  28310 
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1  ARMY  G1 
(PDF)  DAPE  MR  B  KNAPP 

300  ARMY  PENTAGON  RM  2C489 
WASHINGTON  DC  20310-0300 

ABERDEEN  PROVING  GROUND 

32  DIR  USARL 
(20  HC,  RDRLHR 
12  PDF)  L  ALLENDER 

P  FRANASZCZUK 
K  MCDOWELL 
RDRL  HRM 

P  S  A V AGE-KNEPS  HIELD 
RDRL  HRM  AL 
C  PAULILLO 
RDRL  HRM  B 
J  GRYNO VICKI 
RDRL  HRM  C 
L  GARRETT 
RDRL  HRS 
J  LOCKETT 
RDRL  HRS  B 
M  LAFIANDRA 
RDRL  HRS  D 
P  FEDELE  (1  PDF,  20  HC) 

A  SCHARINE 
RDRL  HRS  E 
D  HEADLEY 

30  COMMANDER 
(HC)  US  ARMY  PUBLIC  HEALTH 
CMND 

MCHB  -IP-MHC 
M  ROBINETTE  (20  HC) 

C  JOKEL(IOHC) 

5158  BLACKHAWK  RD 
BLDG  E1570  ROOM  CW-213 
APG  MD  21010 

1  AFRL  RHCB 
(HC)  B  HOBBS 

2610  7TH  ST 
WPAFB  OH  45433 

1  UNIV  OF  CINCINNATI 

(HC)  MECHANICAL  ENGRNG  DEPT 
J  KIM 

2600  CLIFTON  AVE 
CINCINNATI  OH  45221-0072 

1  BOSE  CORP 
(HC)  M  KOSAKOWSKI  MS  271E 
145  PENNSYLVANIA  AVE 
FRAMINGHAM  MA  01701 


1  US  ARMY  MEDICAL  RSRCH  AND 
(HC)  MATL  CMND 
MCMR  RTB  M 
J  LEGGIERI 

FORT  DETRICK  MD  21702-5012 

1  UNIV  OF  ALABAMA 
(HC)  S  A  MCINERNY 

1530  3RD  AVE  S  BEC  356D 
BIRMINGHAM  AL  35294-4461 

1  AFRL/  7 1 1  HPW/RHCB 

(HC)  R  MCKINLEY 

2610  SEVENTH  ST 
WPAFB  OH  45433-7901 

1  NATL  INST  FOR  OCCUPTNL 
(HC)  SFTY  AND  HLTH 

HEARING  LOSS  PREVENTION  TEAM 
W  J  MURPHY 

4676  COLUMBIA  PKWY  MS  C-27 
CINCINNATI  OH  45226-1998 

1  PATTERSON  CONSULTING 
(HC)  J  H  PATTERSON 
11607  BIG  CANOE 
JASPER  GA  30143 

4  AUDITORY  HAZARD  ANLYS 
(HC)  G  R  PRICE 
PO  BOX  368 

CHARLESTOWN  MD  21914 

1  BOSE  CORP 

(HC)  W  RABINOWITZ  MS  27 1 E 
145  PENNSYLVANIA  AVE 
FRAMINGHAM  MA  01701 

1  MCHK  DSH 
(HC)  M  RAMSEY 

1  JARRETT  WHITE  RD 
TRIPLER  AMC  HI  96859 

3  US  ARMY  AEROMEDICAL  RSRCH 
(HC)  LAB 

E  REEVES 
W  AHROON 
K  CASTO 

6901  ANDREWS  AVE 

PO  BOX  620577 

FORT  RUCKER  AL  36362-0577 
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2  L-3  JAYCOR 

(HC)  J  STUHMILLER 
P  CHAN 

10770  WATERIDGE  CIR  STE  200 
SAN  DIEGO  CA  92121 

1  AFRL  RHCB 
(HC)  R  WILDE 

2610  7TH  ST 
WPAFB  OH  45433 

1  JOINT  NON  LETHAL  WEAPONS 
(HC)  DIRECTORATE 
B  A  WRIGHT 
3097  RANGE  RD 
QUANTICO  VA  22134 

1  CENTER  FOR  HEARING  &  DEAFNESS 
(HC)  R  SAL VI 

3435  MAIN  ST 
37  CARY  HALL 
UNIV  AT  BUFFALO 
BUFFALO  NY  14214 

1  CAMPANELLA  ASSOC 
(HC)  A  CAMPENELLA 

3201  RIDGEWOOD  DRIVE 
COLUMBUS  OHIO  43026 

1  WESTERN  MICHIGAN  UNIV 
(HC)  G  A  FLAMME 
ASSOC  PROF 
DEPT  OF  SPEECH 
PATHOLOGY  AND  AUDIOLOGY 
1903  W  MICHIGAN  AVE. 
KALAMAZOO  MI  49008-5355 

1  LUZ  SOCIAL  AND  ENVIRON  ASSOC 
(HC)  G  LUZ 

4910  CROWSON  AVE 
BALTIMORE  MD  21212 

1  ARNOLD  ENG  DEV  CNTR 
(HC)  NAVAIR 
V  BJORN 

740  FOURTH  STREET 
ARNOLD  AFB  TN  37389-6000 

1  N A V AL  S  URF ACE  WARFARE  CTR 

(HC)  DAHLGREN  DIVISION 
W  TRIPP 

18444  FRONTAGE  ROAD  SUITE  321 
DALHGREN  VA  22448-5161 


1  MARCORSYSCOM  MERS 
(HC)  S  TORFIN 

2200  LESTER  STREET 
QUANTICO  VA  22134 

1  ALAN  W  POSTON 
(HC)  8622  PILSEN  ROAD 

RANDALLSTOWN  MD  21133-4432 

1  SEA  05S2 
(HC)  J  MALLINO 

BLDG  197  CUBE  2E-3017 
1333  ISAAC  HULL  AVE  SE 
WASHINGTON  NAVY  YARD  DC 
20376-1210 

1  NAVAL  SAFETY  CTR  LIAISON  OFC 
(HC)  M  GEIGER 

(OPNAV  N09FB) 

2000  NAVY  PENTAGON  ROOM  5C256 
WASHINGTON  DC  20350-2000 

1  PATUXENT  RIVER  NVL  AIR  STA 
(HC)  A  AUBERT 

BLDG  106  RM  225 
22195  ELMER  RD 
PATUXENT  RIVER  MD  20670 

1  OFC  OF  NVL  RSRCH 

(HC)  WARFIGHTER  PERFORMANCE  DEPT 
K  YANKASKAS 
(RM  1033) 

875  NORTH  RANDOLPH  ST 
ARLINGTON  VA  22203-1995 

1  U  S  AFS  AM/OEHT 

(HC)  P  CROWLEY 

2510  FIFTH  ST  SOUTH  BLDG 
W-328 

WRIGHT  PATTERSON  AFB  OHIO  45433 

1  PRGM  MGR  ADVANCED 
(HC)  AMPHIBIOUS  ASSAULT 
S  DROPINSKI 

WORTH  AVE  TECHLGY  CTR 
14041  WORTH  AVE 
WOODBRIDGE  VA  22192 

1  USAFSAM/OEHR 
(HC)  A  T  WELLS 

2510  FIFTH  ST  BLDG  840 
W433.42 

WRIGHT  PATTERSON  AFB  OH 
45433-7913 
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1  NAVAIR  HUMAN  SYSTEMS 

(HC)  ENGRG 
D  GROSS 
MS 5  B2187  R2240 
481 10  SHAW  RD 

PATUXENT  RIVER  MARYLAND  20670 

1  NVL  SURFACE  WARFARE  CTR 
(HC)  DAHLGREN 
A  BERNSTEIN 

G72  COMBAT  SYST  SAFETY  BRANCH 
5375  MARPLE  RD  STE  153 
DAHLGREN  VA  22448-5155 

1  NAVY  MARINE  CORPS  PUBLIC 
(HC)  HEALTH  CTR 
G  F  HOURANI 

620  JOHN  PAUL  JONES  CIRCLE 
STE  1100 

PORTSMOUTH  VA  23708-103 

1  T  LETOWSKI 
(HC)  2015  FAIRWOOD  LANE 

STATE  COLLEGE  PA  16803 

1  MILITARY  OP  MEDICINE  RSCH  PRGM 
(HC)  R  SHOGE 

504  SCOTT  DR 
BLDG  722  RM  43 
FREDERICK  MD  21702 

2  B  AMREIN 

( 1  PDF,  1213  AMBRIDGE  RD 
1  HC)  BEL  AIR  MD  21014 
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1  FRENCH-GERMAN  RSRCH 
(HC)  INSTITUTE  OF  SAINT  LOUIS 
K  BUCK 
PO  BOX  70034 
F-68301  SAINT  LOUIS 
CEDEX  FRANCE 

1  FRENCH-GERMAN  RSRCH  INST  OF 
(HC)  SAINT  LOUIS 
A  DANCER 
PO  BOX  70034 

F-68301  SAINT  LOUIS  CEDEX  FRANCE 

1  CAPTAIN-COMMANDER 
(HC)  F  DE  BOODT 

DEPT  WELL  BEING 
ACOS  WB 

QUEEN  ASTRID  BARRACKS 
BRUYNSTRAAT 
B  1 120  BRUSSELS 
BELGIUM 

1  DEFNC  RSRCH  AND  DEV  CANADA 
(HC)  E  DROLET 

5262  RUE  SNOWDON 
MONTREAL  CANADA  H3W2G1 

1  FRENCH-GERMAN  RSRCH  INST  OF 
(HC)  SAINT-LOUIS 
P  HAMERY 
PO  BOX  70034 

F-68301  SAINT-LOUIS  CEDEX  FRANCE 

1  QINETIQ 
(HC)  S  JAMES 

BLDG  A5  RM  2067 
CODY  TECHNOLOGY  PK 
IVELY  ROAD  FARNBOROUGH 
HANTS  GU14  0LX  UK 

1  CENTRE  DE  RECHERCHES  DU 
(HC)  SERVICE  DE  SANTE  DES  ARMEES 
AGNES  JOB  PH  D 
24  AVENUE  DE  MAQUIS  DU 
GRESIVAUDAN 
F-38702  LA  TRONCHE  FRANCE 

1  P  VAN  DER  VEKEN 
(HC)  ENT  DEPT  SFA  -  CAE 
GEMSTRAAT  61 
9420  ERPE-MERE 
BELGIUM 


20 


1  DIRECTORATE  OF  FORCE  HEALTH 
(HC)  PROTECTION 

CANADIAN  FORCES  HEALTH 
SERVICES  GROUP  HEADQUARTERS 
NATIONAL  DEFENCE 
S  TSEKREKOS 
1745  ALTA  VISTA  DRIVE 
OTTAWA  ONTARIO  CANADA 
K1A  0K6 

1  AMNON  DUVDEVANY  PHD 
(HC)  HEAD  OF  ACOUSTICS 
I  D  F  MEDICAL  CORPS 
MILITARY  PO  BOX  02149 
I  D  F  ISRAEL 
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